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The relation between the boolean functions and Bell inequalities for qubits is analyzed. The 
connection between the maximal quantum violation of a Bell inequality and the nonlinearity of 
the corresponding boolean function is discussed. A visualization scheme of boolean functions is 
proposed. An attempt to classify Bell inequalities for qubits is made, a weaker result (classification 
with respect to Jevons group) is obtained. The fractal structure of the classification is shown. All 
constructs are illustrated by Mathematica code. 
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I. INTRODUCTION 

In their famous paper [l[ Einstein, Podolsky and Rosen 
(EPR) suggested a Gedankenexperiment which, as they 
believed, must prove the incompleteness of quantum me- 
chanics. An interesting analysis of this problem wasgiven 
by Bohr 0. Some progress was achieved by Bell [3[. He 
showed that under assumption of the EPR arguments 
some inequalities must be fulfilled. If Bohr's arguments 
are correct then these inequalities can be violated. It 
was only in 1980 when the Bohr's arguments were exper- 
imentally verified 0, H, [E 0, @] • Now the arguments by 
EPR are considered to be incorrect, but nevertheless it 
is the work [l[ that initiated the discussion of basics of 
quantum mechanics. 

In this work I analyze the connection between the 
boolean functions theory and Bell inequalities for multi- 
qubit systems (which were obtained in Q). Surprisingly 
enough, in many aspects Bell inequalities theory is anal- 
ogous to the combinatorial problems of computer logic 
circuits design developed a year earlier then the Bell's 
work [|| appeared. There is also a relation between Bell 
inequalities and applications of boolean functions theory 
to classical cryptography. For example, classification of 



is closely connected to 
|, [ll[ , and the maximal 



Bell inequalities discussed in 
the Jevons group, studied in 
quantum violation of a given Bell inequality is connected 
to the nonlinearity of the corresponding boolean func- 
tion. I made an attempt to clarify these connections. But 
there are still many open questions, both combinatorial 
(like classification of Bell inequalities with respect to the 
group Q n ) and analytical (like calculation of the maxi- 
mal quantum violation Vf). Another interesting problem 
is the connection between the maximal quantum viola- 
tion and the uncertainty relation for boolean functions. 
In this work it is shown that the Bell inequalities whose 
maximal quantum violation is the largest (Mermin in- 
equalities), minimize this uncertainty relation. 

All mathematical constructions discussed in the text 
are illustrated with Mathematica. I choose Mathemat- 
ica since it has an extremely flexible and unified pro- 
gramming language and a very rich set of built-in math- 
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cmatical functions. Due to this it is possible to present 
the algorithms illustrating the discussed quantities in a 
very compact form. Using Mathematica, it is possible to 
code all the illustrating examples using high-level con- 
structions, avoiding worrying about low-level program- 
ming details which have nothing to do with the problem 
under study. I do not pretend to give the most effec- 
tive Mathematica code for calculating different features 
of boolean functions, my goal is to present a compact 
and ready-to-use working code. All the examples can be 
typed in and run provided that they are entered to Math- 
ematica in the given order. I also presented the C source 
code for the fast Walsh-Hadamard transform and a way 
to turn it into an executable program which can be used 
in Mathematica. 



II. BOOLEAN FUNCTIONS 

The Bell inequalities for a multi-qubit system, which 
were obtained in are closely connected with boolean 
functions theory. There is a natural one-to-one corre- 
spondence between the set of all B n = 2 2 boolean func- 
tions of n boolean variables and the set of Bell inequali- 
ties for n-qubits. In this section I give a short overview 
of the notions and results of the boolean functions the- 
ory which are needed for applications to Bell inequalities. 
The most important notion discussed in this section is the 
Walsh-Hadamar transform, which turns out to be the co- 
efficients of the Bell inequality corresponding to a given 
boolean function. 

Here I also introduce a visualization technique of 
boolean functions which is useful to graphically represent 
different classes of boolean functions. This approach is 
based on the fact that for all n the number of boolean 
functions of n variables is a square (in fact, B n = B\_^), 
so that one can associate the boolean functions with the 
cells of a B n -i xB n _i array. Briefly speaking, this is done 
numerating boolean functions with integers, interpreting 
the vectors of their values as binary decompositions and 
then using division modulo B n _\. In this section I visu- 
alize boolean functions with respect to their degree and 
uncertainty. In both cases the pictures show some kind 
of fractal behavior. 
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A. Boolean vectors 

Let Z2 = Z/2Z be the finite field with two elements 
Z 2 = {0, 1}. The sum and the product of elements a, b <G 
Z2 are denoted as a © b and ab respectively. The product 
of n copies of Z2 we denote as V n = ZJ? and refer to 
its elements as (n-dimcnsional) boolean vectors. It is 
clear that \V n \ = 2". The notation (x, y) is used for the 
scalar product of two boolean vectors x = (xi, . . . , x n ) 
and y = . . . ,y n ): 

(x,y) =xij/i © ... ®x n y n . (1) 

There is a natural one-to-one correspondence b n be- 
tween the set V n and the set B n = {0, . . . , 2™ — 1}: 

n 

b n :V n 3x = (x u ...,x n ) -> x = ^2x i 2 n ~ i £ B n . (2) 

i=l 

In other words, the boolean vector x corresponds to the 
integer x whose binary representation is given by x. The 
most significant bit in the binary decomposition of x is 
the first component of x and the least significant bit is 
the last component. The correspondence b n is natural 
in the following sense. For an integer m < n the set 
V m can be identified with a subset of V n by padding Tri- 
dimensional boolean vectors on the left to extend them 
to n dimensions. Then the diagram 



n n (3) 

is commutative. This means that one can apply b n to 
m-dimensional boolean vectors with m < n. The corre- 
spondence b n is illustrated by Table HI 



X 


6n(x) 


m = 1 


m = 2 


m = 3 


(0,0,. ..,0,0,0) 





61 (0) 


62(0,0) 


63(0,0,0) 


(0,0,. ..,0,0,1) 


1 


Mi) 


62(0,1) 


63(0,0,1) 


(0,0,. ..,0,1,0) 


2 




62(1,0) 


63(0,1,0) 


(0,0,. ..,0,1,1) 


3 


62(1,1) 


63(0,1,1) 


(0,0,..., 1,0,0) 


4 




63(1,0,0) 








(0,0,. ..,1,1,1) 


7 


63(1,1,1) 








(1,0,. ..,0,0,0) 


2 n-l 






(1,1,..., 1,1,1) 


T - 1 



TABLE I: The correspondence b n between boolean vectors 
x £ V n and integer numbers x £ B n . 

The set V n can be ordered in many ways. I use the 
lexicographical order: for x = {x\, . . . , x n ) and y = 
(yi, ■ ■ ■ , Vn) the notation x < lcx y means that there is 



k, 1 ^ k ^ n, such that xi = yi for i = 1, . . . , k — 1, but 
Xk = and yk = 1- Here and 1 are elements of Z 2 , 
but if one considers them to be integers, the same can 
be shorter formulated as: x < lcx y if and only if the first 
non-zero difference y\ — X\, . . . , y n — x n is positive. Since 

fc-i 

^2 i = 2 fc -K2 fc , (4) 

i=l 

it is clear that x < Icx y if and only if a; < y, where 
x = b n (x) and y = b n (y). This means that b n preserves 
the order: 

(K,< lcx )^ (£„,<). (5) 

There are two very important functions on V n : Hem- 
ming weight wt(x) and Hemming distance d(x, y). The 
Hemming weight wt(x) of a boolean vector x G V n is 
defined to be the number of non-zero components of 
x. The Hemming distance between two boolean vec- 
tors x, y £ V n is the number of positions where com- 
ponents of x and y differ. It is clear that the distance 
can be expressed in terms of the Hemming weight as 
d(x, y) = wt(x©y). Both the notions play an important 
role in different applications of boolean functions theory. 

B. Boolean functions 

A boolean function /(x) = f(x\, . . . , x n ) of n boolean 
variables is a map / : V n — ► Z2. The set of all boolean 
functions of n boolean variables we denote as F n . There 
is a one-to-one correspondence between F n and V^n : 

f»3/- (/(6n 1 (2 B - f(KHo))) e V*». (6) 

Due to this we have ji^l = 2 2 . 

Any boolean function / € T n can be represented in 
the following form: 

/(x) = g f (yW, (7) 
yev„ 

where gf G F n and x y = xf 1 . . . x n n (with 0° = 1) or in 
the equivalent form as 

n 

/(x)=c o ©0 c il ...i h x il ...x ik . (8) 

k—l l^ii <.. .<2fc 

The relation |7]) can be inverted and the function g/ is 
expressed through / as 

5/(y) = 0/(x), (9) 

where x =<: y means that x^ ^ for k = 1, . . . , n (note 
that =<! is not a linear order). From this it follows that 
the relation / <-> gf is one-to-one, which is referred to 
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((-l) <U ' V> )u.v 



(16) 



as Mobius transform. It is an involution: if gt = h then The tensor product H n = H® n of n copies of H reads as 
9h = f, or g gf = f for all / G F n . 

The degree dcg / of a boolean function / G F n is de- 
fined to be the maximal number d such that there is 
y G V n with wt(y) = d and g/(y) = 1. Boolean func- 
tions of degree 1 are called affine; they can be represented 
as follows: 



For any / G F n we define two boolean vectors W/, z/ G 
via 



a(x) = cq © CiXi 



c n x n = c © (c,x), 



(10) 



where c = .9/(0, . . . , 0), = 5/(0, . . . , 1, . . . , 0) (1 is on 
the i-th position) and c = (ci, . . . , c n ). The set of all 
affine functions is denoted as A n . If Co = then the affine 
function is called linear; such functions are denoted as I 
(so that Z(x) = (c,x)). The set of all linear functions is 
denoted as L n . It is clear that \L n \ = 2 n and \A n \ = 
2 n+1 . 

A boolean function / G F n is called homogeneous of 
degree k if <7/(y) = for all y G V n with wt(y) ^ k. If, 
in addition, <?/(y) = 1 for all y G V n with wt(y) = k, 
then / is referred to as the fc-th symmetric function and 
denoted as Sk- 



Sfc(xi, . . . ,x n ) 



e ■ 

l^il<...<ifc^n 



(11) 



For boolean functions one can introduce the Hemming 
weight and distance in the same way as it was done for 
boolean vectors: wt(/) is defined to be the number of 
x G V n with /(x) = 1, and d(f,g) is the number of 
x G V n with /(x) 7^ ff(x). A boolean function is called 
balanced if wt(/) = 2"~ 1 . If we denote / © g G F n the 
point-wise sum of / and g, 



(/©ff)(x) = /(x)© 5 (x), 



(12) 



then d(f, g) = wt(/ffip). The distance d(f, M) between a 
boolean function / G F n and a nonempty subset M C J 7 ,, 
is defined via 



d(/,M) = min </,<?). 



(13) 



The distance TV/ between / and A n is called nonlincarity 
of /: Nf = d(f, A n ). It is a very important characteristic 
of boolean functions. 



C. Walsch-Hadamar transform 

In studying properties of boolean functions the no- 
tion of Walsch-Hadamard transform can be very useful. 
The Walsch-Hadamard transform of a boolean function 
/ G F n is the integer-valued function Wf of n boolean 
arguments defined via 



W f (u) = £ (-1) 



/(x)e<x,u> 



Let us introduce the Hadamard matrix H, 



H 




(14) 



(15) 



w / = (W>(u)) u6Vn> z f = ((-1)/W) xev „. (17) 

The Walsch-Hadamard transform (|14[) can be written in 
the following compact form: 



H„z 



nZf. 



(18) 



To calculate any component of w/ according to the def- 
inition (fl"4"|) it takes 2™ — 1 additions; to calculate all 
2™ components it takes 2™(2" - 1) - 2 2 ™additions. Us- 
ing special properties of the matrix H n , from (|18p it is 
possible to calculate the vector w/ (given z/) with only 
n2 n additions (fast Walsch-Hadamard transform). The 
algorithm and its realization in C (to be used in Mathe- 
matics) are described in Appendix B. 

Note that all components of w/ are even, since a sum 
of an even number of equal terms is even. As an example, 
let us calculate w; for a linear function Z(x) = (c, x). We 
need the following relation: 



V (-1)<*.»> = 2"5 X , . 
uev n 



(19) 



This relation is obvious due to the equality 

71 

- (_!)(«,«) = JJ(i + (20) 



u<EV„ 



fc=l 



Using the relation (|19p from the definition (jT4j) one can 
easily get 



W,(u) = 2 n 6 u , c , w / =2"e c 



(21) 



where e c = (0, . . . , 1, . . . , 0) (1 is on the 6„(c)-th posi- 
tion). 

The Walch-Hadamard transform (|T4|l is invcrtible; the 
inverse transform reads as 

(-!) /(x) 4^ (-l) <X ' U> ^/(u), (22) 



or, in matrix notation, as 

z f = 2~ n H n w f . 



(23) 



This relation is a trivial consequence of (fT5)) due to the 
simple fact that H% = 2 n E 2 ^, or H~ x = 2~ n H n . Using 
the same algorithm of fast Walch-Hadamard transform, 
one can quickly calculate z / given w / . 

Below we will need the following statement: for an 
integer-valued function W(u), u G V n there is f G F n 
such that Wf(u) = ^(u)^ if and only if W(u) satisfies 
the condition 



J2 W(u)W{u®v) 



uev„ 



if v = 0, 
if v^O. 



(24) 
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In particular, the Walsch-Hadamar transform Wf(u) of 
any boolean function f G F n satisfies the Parseval equal- 
ity: 



Wf(u) = 2 2n . 



(25) 



uev„ 



The proof of this statement is quite simple. The Walsh- 
Hadamard transform of any boolean function / G J- n 
satisfies the condition (f2~4"|) , which can be easily derived 
from the definition (fl4|) using the equality (fT9|) . The less 
trivial part is to prove that the condition (|24|) is sufficient. 
It is enough to prove that the function 



n*) = 4 E (-i) <x>u> ^(u) 



(26) 



takes only two values ±1. In other words, it is sufficient 
to prove that F 2 (x) = 1 for all x G V n . We have 

p2 ( x ) = ^ E (-l) (x ' u>e(x ' v> ^(u)^(v). (27) 
u,vev„ 

For any fixed u G V n the inner summation over v G V n 
is equivalent to the summation over u © v: 



f2 W = ^t E (-i) (x ' v> ^(u)^(uev). (28) 

u,vEV„ 

Changing the summation order we get 

=iE (-!) <X ' V> E W(u)W(u®v). (29) 



v£V„ 



u£V n 



Due to our assumption the inner sum differs from zero 
only for v = and finally we get F 2 (x) = 1. This 
completes the proof. 

For nonlinearity one can get the following result: 



JV/ = 2"- ] --mnxjW f (u)\. 



1 



(30) 



Let us denote NW f = \{u e V n \W f {u) ^ 0}\, the num- 
ber of non-zero components of w/. Due to the Parseval 
equality (|25|) we have 



(31) 



2 2 " = E w H u ) < NW f™% \ w f( u )\ 



and from (|30p we get the inequality 

2 «-i 



N f < 2 r 



(32) 



/ 



For example, for a linear function I G L„ from (|21|) we 
have AW/ = 1 and from ([32]) it follows A ; = 0. 

It is also possible to get an absolute upper bound for 
nonlinearity. According to the Parseval equality (f25j) we 
have 



max|W / (u)|^2"/ 2 , 

UE V n 



(33) 



and due to the relation (|30"|) we get 



N f sC 2' 1_1 - 2™/ 2_1 . 



(34) 



As we will see, for an even n this bound is exact. For an 
odd n the exact bound is unknown. 



D. Walsh-Hadamard transform as a representation 

Now another point of view on the Walsch-Hadamard 
transform will be presented. From the definition (jT4j) one 
can easily derive the following equality: 

J2 VF/(uew)W 9 (wev) = 2' l W7e 9 (u®v), (35) 

wEV„ 

valid for all f,g G F n . In partial case of / = g it reads as 
W)(uew)W/(wev) = 2 2n ^ v . (36) 

wev„ 

Let us for any / G F n introduce the 2™ x 2"-matrix Wt 
via 



Wf = ^(^/(Uffiv)) u , v£ y„ 



(37) 



The equality (|36|) means that any matrix Wf is a nontriv- 
ial root of the identity matrix: W 2 = E 2 ™ ■ The equality 
(|3"5"j) means that for all /, g € F„ 



W/W fl = W7©„- (38) 

In other words, the map IF, 

VF:F„^GL(R 2 "), f^Wf, (39) 

is a representation of the additive group F n in R 2 . Its 
character is given by 



X (/)EETr(IF / ) = (-ir t ^. 



(40) 



Since F n is a commutative group, W is equivalent to 
a direct sum of one-dimensional representations. This 
decomposition can be given explicitly First, let us find 
eigenvalues an eigenvectors of Wf. Namely, let us show 
that for all w G V n the & n (w)-th column z w of H„ , 



f)=((-l) M )uEK, 



(41) 



is an eigenvector of Wf with the corresponding eigenvalue 

(_l)/(w). 



(wy*(w)) u = J2 w>(uev)(-i)< v - 

v£V n 

= E (-i) (u0v - w> ^/(v) 



(42) 



v€V„ 



(-l)JW(-l) 



(-lK< w >z(w) u . 
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11) 



H 



h) 



If ffi / (x) 



jj®n — |W/) 



FIG. 1: Figure 

This means that the matrices Wf are diagonal in the 
basis {z(w)} we v„. In particular, for the determinant of 
Wf we have 

detWf = (-l) wt{f) . (43) 
The diagonalization of Wf reads as 

W f = 2~ n H n diag((-l)£WjfT n . (44) 

For any fixed w € V n the map 

/ - (-l) /(w) (45) 

is a one-dimensional representation of -F n and the equal- 
ity (|4"4"| gives the decomposition of into the direct sum 
of such representations. The line spanned by z(w) is an 
invariant subspace of R 2 on which W acts as ([45)) . 



E. Autocorrelation function 

Another important notion is the autocorrelation func- 
tion A/ iS (u) of two boolean functions f,g € F n , which is 
defined as follows: 



/(x)es(x©u) 



xey„ 



± (-l)^W f (v)W g (v). 



(46) 



v£V„ 



Since the map x — > x © u is one-to-one for a fixed u and 
idempotent, i.e. (x u) © u = x for all x and u, the 
autocorrelation A/ j9 (u) is symmetric with respect to / 
and g: A/. s (u) = A g j(u). For f = g the autocorrelation 
A/j(u) = A/(u) is referred to as the autocorrelation of 
the function /. 



F. Physical interpretation 

Here another interpretation of Walsch-Hadamard 
transform is presented. Consider an n-qubit system. For 
any / £ F n let us introduce the state \ Wf) via 



\ W f) = 4 £ W f (u)\u). 



(47) 



The state \Wt) can be obtained using the circuit shown 
in Fig. [T] Here the Hadamard gate H is given by the 



matrix H = 2 X / 2 H and the in + l)-gubit gate Uf acts 



as 



C/ / (|x)|y)) = |x)|y©/(x)). 



(48) 



Straightforward algebraical manipulations show that if 
the input state is |0)®"|1) then the output state is 
IW 7 /)!— ), where 



|±> = -^(|0>±|1». 



(49) 



Due to the Parseval equality (|2"5j) the states |W/) are 
normalized: 



(50) 



ugV„ 



but they are not orthogonal: 



From this we can conclude: \Wf) and | W g ) are orthogonal 
if and only if A/ iS (0) = 0. 

For any single-qubit operator A and for any v = 
(vi, . . . , v n ) S 14, we use the notation 



A v = A 1 



A v - 



(52) 



for the n-qubit factorizable gate, with A:-th component to 
be A (if Vk = 1) or 1 (if v*, = 0). The operator with the 
matrix 



X 



1' 

1 o, 



(53) 



is called NOT gate; it acts as X\0) = X\l) = |0). It 
is clear that for any v G V n 

A v |u) = |uffiv). (54) 

From the equality we have 

(W f \X v \Wf) =^J2 W7(u)W>(uffiv) = S v , , (55) 
uey„ 

i.e. all diagonal matrix elements of the from (W r /|A V |W/) 
of any non-trivial operator (with v^O) operator X v are 
equal to zero. 

Let us also consider the phase gate 



Z 



1 
.0 -1 



(56) 



It acts as Z\0) = |0), Z\l) = It is clear that 

Z v |u) = (-l)^|u). (57) 
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There is a generalization of the equality ([51]) : 

(W f \Z v \W g ) = ^A /iS (v), (58) 

valid for all v G V n . 

For any state \Wf) let us consider 2™ states 

X v \W f ) = ±J2W f (u®v)\u). (59) 

Below it will be shown that 

X v \W f ) = \W U ), (60) 

where / v (x) = /(x) © (x, v). Due to the relation (fBTj) 
the states |W/ V ) and |W/ W ) with different v and w are 
orthogonal: 

(W fv \W f J = 6 v , vr . (61) 
One can easily prove that for all v G V n 

E W f (u®v) = E W M) = (-l) /(0) 2", (62) 
from which we get the following relation 

E iw/v) = (-i)^ ' E i u >- ( 63 ) 

v6V„ uev„ 

Since the matrix Wf is non-degenerate according to (|43| , 
for any fixed / G ^ the 2" vectors |W/ V ), v£V„ form 
a basis of the state space of n qubits. The equality ([63]) 
shows that the sum of the new basic vectors (the diagonal 
of the parallelepiped spanned by the new basic vectors) 
is modulo the sign equal to the old one. 

G. Uncertainty relation 

Let NAf be the number of u G V„ with A/(u) ^ 0. 
Then the following statement is valid : {or all f G F n 
the numbers NWf and NA f satisfy the inequality 

NWjNAf > 2™. (64) 

The quantity U(f) 1, defined via 

U(f) = ^NWfNA f , (65) 

is referred to as the uncertainty of / G F n . 

In cryptographical applications of boolean functions 
theory the numbers NWf and NAf play an important 
role. In some applications it is necessary to use boolean 
functions / with small NWf , in others with small NA / . 
The inequality (f64|) shows that both the numbers NWf 
and NAf cannot be small, they are subject to jSJ). In 
this sense the inequality can be called the uncer- 
tainty relation for boolean functions. 



f(xi,x 2 ) 


(0,0) 


(0,1) 


(1,0) 


(1,1) 




(hi) 

















1 


(1,1) 


(l + a;i)(l + X2) 


1 











2 


(1,2) 


X2 + X\X2 





1 








3 


(1,3) 


X2 


1 


1 








4 


(1,4) 


XI + X\X2 








1 





5 


(2,1) 


X\ 


1 





1 





6 


(2,2) 


X\ + X 2 





1 


1 





7 


(2,3) 


1 + X1X2 


1 


1 


1 





8 


(2,4) 


X1X2 











1 


9 


(3,1) 


I+X1+X2 


1 








1 


10 


(3,2) 


l + Xi 





1 





1 


11 


(3,3) 


1 + XI + X\X2 


1 


1 





1 


12 


(3,4) 


1+X 2 








1 


1 


13 


(4,1) 


1 + X2 + X\X2 


1 





1 


1 


14 


(4,2) 


X\ + X2 + X\X2 





1 


1 


1 


15 


(4,3) 


1 


1 


1 


1 


1 


16 


(4,4) 



TABLE II: Boolean function visualization for n = 2. 



H. Visualization of boolean functions 

Now an approach to the visualization of the set of 
boolean functions will be presented. The set F n can be 
identified with the set #2™ using the map b n (remember 
© and ©): 

B n = b 2 n : F n 3 f E /(x)2 Mx) G B 2 n. (66) 

x£V„ 

Then with each function / G F n one can associate a pair 
of integers (i,j) such that 

i?„(/) = (*-l)2 2 "" + (j-l). (67) 

When / varies over T n , the corresponding pair (i,j) runs 
over the square 

Sq„ = {(z,j)|l^*,J^2 2 ""}. (68) 

As an example, the correspondence between / G F 2 and 
Sq 2 is shown in the Table HT1 This square can be depicted 
as a 2 2 " 1 x 2 2 " 1 array of cells, with each cell marked 
with a definite color. For example, one can mark the 
cells corresponding to the functions of the same degree 
with the same color. Below we present pictures of F n col- 
ored with respect to different characteristics of boolean 
functions, not only degree. 

Let us start with the visualization of boolean functions 
of different orders. In Fig. [2] the cases of n = 2, n = 3 
and n = 4 are shown. Each subfigure contains n + 1 
different colors. Then let us visualize boolean functions 
in the three cases with respect to the uncertainty U(f), 
(frJ5"j) . Fig. [3] shows the same three cases of n — 2, n = 3 
and n = 4. 
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(a)n = 2. White corresponds to U(f) = 1. (b)n = 3. Light gray corresponds to U(f) = 1, 

black to 8. 




(c)n = 4. White corresponds to U(f) = 1, yellow to 35/8, red to 8 and blue to 16. 
FIG. 3: Uncertainties for three different cases. 



9 



III. BELL INEQUALITIES 

In this section the Bell inequalities for n qubits and 
especially the extremal subclass of Mermin inequalities 
are considered and the relation between the maximal vi- 
olation of a given Bell inequality and the nonlincarity of 
the corresponding boolean function is discussed. 



A. Construction of Bell inequalities 

Consider n pairs of random variables Ai(0), Ai(l), 
i = l,...,n taking only two values ±1. Let E(u), 
u = (u\, . . . ,u n ) G V n be the mathematical expectation 
of the product of n variables Ai(ui), i = 1, . . . , n: 



E(u)=M(A 1 (u 1 )...A n (u n )). 



(69) 



Clearly, all 2™ expectations E(u), u G V n , satisfy the 
inequality 



\E(u)\ < 1. 



(70) 



Can any 2™ numbers E(u) subject to ((70|) be the math- 
ematical expectations according to (|6"9"|) ? The answer is 
negative. 

Let us illustrate this statement by a simple example in 
the case of n = 2. For E{u\,u 2 ) we have 

E( Ul ,u 2 ) = 2P(Ai(u x ) - A 2 (u 2 j) - 1. (71) 

From this expression one can conclude that 

E{u u u 2 ) = 1 & P(Ai(«i) = A 2 {u 2 )) = 1, (72) 

and that 

E( Ul , « a ) = -1 P(Ai(«i) = A 2 (u 2 )) = 0. (73) 

It is easy to see that the numbers (1,1,1,-1) cannot 
be the expectations (E(0,0), E(0,1), E(1,0), E(l,l)). In 
fact, from fl72J) for (ui,u 2 ) = (0,0), (0,1) and (1,0) it 
follows that in such a case for any realization of Ai(Q), 
A 2 {0) and A 2 (l) we would have Ai(l) = A 2 {0) = 
AiiO) = A 2 (l), but from (JT3J) for (ui,ti 2 ) = (1,1) we 
would have A\(l) ^ A 2 (l). This contradiction proves 
that the numbers (1,1,1,-1) cannot be the expectations 

(59). 

Now a set of necessary and sufficient conditions for 
E(u), u G 14 to be the expectations (|6"9"|) is derived. 
These conditions were obtained in @ . Our approach fol- 
lows the idea of [l4|. For a fixed / G F n consider the 
random variable Af defined via 



.4 



(-1)^11(^(0) + (-1)^(1)). (74) 



i=l 



For any realization of Ai{v,j) the product in (|7^|) differs 
from zero only for one x G y n , and for this x the product 



is equal to ±2 ra ; for all other 2 n — 1 boolean vectors x G V n 
it equals to zero. We see that for any / G F n only one 
term in the sum ([7^)1 differs from zero and equals to ±2™ 
so that the sum takes only two values ±2™. From this 
one can conclude that 



\M(A f )\ < T 



(75) 



The random variable Af can be written in the following 
form: 



A f= E(-!) /(x) EC-i^AiM.-.A^) 

xev„ uev„ 
= Yj Wf(u)A 1 (u 1 )...A n (u n ). 



uEV n 



(76) 



Taking the mathematical expectation, the inequality ([75 
becomes 



J2 W f (u)E(u) 



u6V„ 



(77) 



Note that for the function /(x) = 1 © /(x) we have 
Wj(u) = — Wf(u), and due to this the inequality (f77|) 
for f € F n is equivalent to two inequalities 



£ W>(u)S(u)<2 f 



(78) 



for / and /. The 2 2 " inequalities (JTSJ) for / G F n are 
referred to as Bell inequalities. They form a necessary 
and sufficient condition for E(u) to be the mathematical 
expectations (|6"9"|) . 

Note that the inequality ((70|) is equivalent to two in- 
equalities 



E(u) and - E(u) < 1. 



(79) 

The first inequality is the Bell inequality ([78]) with the 
linear function /(x) = Z(x) = (u, x); the second one 
corresponds to the affine function 7(x) = 10 (u, x) , so the 
trivial inequalities (|70[) are Bell inequalities with affine 
functions. 

The Bell inequalities (f78|) can be written in the follow- 
ing formal form: 

(W f \E) < 1, (80) 
where \E) is a (non-normalized) state defined via 

(81) 



\E) = J2 ^( u )i u ) 

u£V„ 



If Ai(ui), i = 1 
then we have 



, n are independent for all itj G 7i 2 , 
^(ui, . ..,u n ) = qi(ui) . . .q n (u n ), (82) 
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where qi(u) = 2P(Ai(u) = 1) — 1, and \E) is factorizablc 

n 

|S)=®(?i(0)|0) + ft (l)|l)). (83) 

i=l 

Since |<?i(u)l ^ 1 it is easy to check that any factorizable 
state (|83|) satisfy all the inequalities ([80]) . But even if 
there are correlations between Ai(v,i), their mathematical 
expectations ([6"9")l satisfy the Bell inequalities 1(75]). 



It was shown that the largest maximal violation reads 

as 

maxjjf = 2 { - n - 1 ^ 2 . (87) 

The inequalities on which this upper bound is attained 
are called Mermin inequalities [131 ]. In the next sub- 
section the boolean functions corresponding to these in- 
equalities are explicitly constructed. 



B. Violations of Bell inequalities 



C. Mermin inequalities 



Now let Ai(ui) be Hermitian operators (observables) 
with the spectra {—1, 1} and for a given n-qubit state g 
let E(u) be the quantum-mechanical average 

E(u) = (A 1 (ui)...A n (u n )} s = Tr(gA 1 (u 1 )...A n (u n )). 

(84) 

If we assume that for the state g the result of the measure- 
ment of the observable Ai(ui) can be described by a ±1- 
valued random variable Ai(ui), then the quantum me- 
chanical average ([84]) equals to the mathematical expec- 
tation l|6"9")) and the quantities E(u) defined in l|84p satisfy 
the Bell inequalities (|78|) . Such a state g is called classi- 
cally correlated. Surprisingly, there are quantum states 
such that the quantities (|84|) violate the Bell inequalities 
(for a definite choice of observables Ai(ui)), which means 
that the correlations of such states are stronger than any 
classical ones. 

It is interesting to find out up to what extent a given 
Bell inequality can be violated in quantum case. In Q 
it was shown that the maximal violation vj of the Bell 
inequality ([78]) corresponding to the boolean function / e 
Fr, reads as 



1 

Vf = — max 

2™ tp lt ...,(p n e[0,2ir] 



E 



W f (u) 



(85) 



where <p = {tpx, . . . ,tp n ) and (<p,u) = J2k=i u kVk- Us- 
ing the definition (fl"4|) of Wf (u) , this expression can be 
rewritten in the following form: 



vr — max 



E 

xev„ 



H) /W IIH\W 



k=l 



(86) 



where t (ip) = cos ip and t\{ip) = snap. In general, it is 
not an easy optimization problem and the explicit form 
of Vf is unknown. The numerically calculated values of 
Vf for the small n are shown in Fig. [4] Studying the 
numerically obtained results one can conclude that there 
is no unique relation between the maximal quantum vi- 
olation Vf and the nonlinearity Nf, but nevertheless the 
higher Nf the larger Vf (at least for n = 2,3,4). In 
other words, the higher nonlinearity of a boolean func- 
tion function the stronger maximal quantum violation of 
the corresponding Bell inequality. 



Let Ai(u), i = 1, ...,71, u G Z2 be ±l-valued ran- 
dom variables. For x,y € Z2 define the random variable 
M n (x, y) via 

n 

M n (x, y) = Im IJ(j4 fc (aO + iA k (y)). (88) 

k=l 

For an odd n the Mermin inequality reads as 

M(M n (0,l)) s=:2("- 1 )/ 2 , (89) 
and for an even n it reads as 

M(M n ^(0,l)(A n (0) + A n {l)) 



+M„_ 1 (1,0)(A„(0)-A„(1))) < 2"/ 2 . 



(90) 



After multiplying by a proper number, these inequalities 
can be written in the form (|78| with Wf(u) = M(u) 
given by 

, , [O if wt(u) is even, . . 

M u) = <^ / , „ +1 91 

1(_1)-4 L 2^- if wt(u) is odd, 

in the case of odd n, and by 

M(u) = ( M(U ' )+M(n ? 02) 
V 1 \M(u')-M(u') u n = l, y ' 

in the case of even n, where u' = (v,\, . . . , u n _i) so that 
M(u') is defined by (f9T]) . 

The function M(u) defined by (|9ip or (|92p satisfies 
the condition (|24[) . Let us first consider the case of odd 
n. For v = we have 

A/2 ( u ) = E 2 " +1 = 2" +1 2"- 1 = 2 2 ". (93) 

u£V„ wt(u) is odd 

For v 7^ first consider the case of odd wt(v). Then 
wt(u © v) is odd for any u with odd wt(u) and due to 
this all terms in the sum in (|24p are zero and the condition 
([24]) is satisfied. The case of even wt(v) is more difficult. 
If wt(u) = 2m + 1 then M(u) ~ (-1)" 1 (we omit the 
common factor 2^ n+1 ^ 2 ) and 



wt(u © v) = wt(u) + wt(v) - 2k, 



(94) 
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where k = wt(u&v). We have 

M(uffiv) ~ (_l)«+wt(v)/2-fc ) ( Q5 ) 

and the sum in the condition reads as 

ueV„ k wt(u&v)=fc 



The number of terms in the internal sum is 

wt ( v )^ 2 n-wt(v)-l 

k 



(97) 



According to our assumptions n is odd and wt(v) is even, 
hence n — wt(v) — 1^0 and the sum can be calcu- 
lated as 

J2 J2 (-l) fc = 2"- wt ( v )- 1 (l - l) wt(v ' = 0. (98) 

k wt(u&v)— k 



This proves that the condition ([24]) is satisfied also in the 
case of even wt(v). The case of odd n was completely 
considered. The proof for the case of even n easily follows 
from thew definition (|9"2"j) . 

Now we will find the boolean functions which corre- 
spond to the Mermin inequalities or, in other words, 
which maximize Vf. First, we will find the boolean 
function m G F n whose Walsch-Hadamard transform is 
M(u). We start with the case of odd n. According to 
we have 



wt(u)-l 



wt(u) is odd 

= 2-^ J2 (-l)^fc((-l) x V ••,(-!)*"), 

k is odd 



(99) 



where Sk is the k-th symmetric polynomial 

8k(zi, . . . ,z n ) = ^2 z n ...z ik . (100) 

Note that the last sum in (l99l) can be written as 



J2 (-l)^s k ( Zl , z n ) = (-1)^ Ref]^ - 

(101) 



l^k^n 
k is odd 



1=1 



and (|99p can be simplified 

n 

(-l) m(x) = -(-2)-TBe JJ(j _(_!)**). (102) 

fc=i 

Using the relation 

<««> 



we get the equality 

(.jjmW = (_!)wt(x)-2=l Rc _ .^i(wt(x)-2 i ^)f^ _ 

(104) 

Explicitly the function m(x) can be written as 



m(x) 



if wt(x) - 2=1 = 0, 3 (mod 4), 

1 if wt (x 



(105) 

2=1 = 1,2 (mod 4). V ' 



For any boolean vector x = (x±, . . . , x n ) € V n we have 
0x fc =wt(x), 0x 4 s j = f W ^ X) ) (mod2), (106) 



k=l i<j x 

from which we get the following result: 



m(x) = < 





B © 


X r^X J 


if n = 


1 (mod 


k=l 










© XiX. 






if n = 


3 (mod 


i<j 

Tl 










1© © 

fe=l 




i<j 


if n = 


5 (mod 


10© 






if n 


7 (mod 



(107) 

This gives the decomposition of m in the case of odd n. 
hi the case of even n from the definition (l9"2"l) we have 



(108) 



(_!)«»(*) =2 -n ^ (-l)( x »(M(u) +M(u)) 

+ 2~ n ^ (-l) (x '' u)+x "(A/(u)-M(u)), 
uev n _i 



Since n — 1 is odd, from the previous case the we can 
conclude that the following relations arc valid: 

J2 (-l) (x '' u) A/(u) = 2"- 1 (-l) m «>, 
uev„_i 

( — l)' X ' ,U 'M(u) = 2"~ 1 (— l) m ( x ')®wt(x)®a;„^ 

uGV„_i 



(109) 



(110) 



from which we get 

(_l) m ( x ) — i(_iy™( x ') 

(l + (_i)»t( x )+^ + _ (-i) wt ( x )), 

what is equivalent to the relation 

. . I to(x') if wt(x) is odd, 

m(x) = < 

ffl(x )®i„ if wt(x) is even. 
Explicitly m(x) reads as 

\ 

m(x) = m(x') © wt(x')x n = m(x') I irfc ] £„. 

(112) 



(111) 
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This gives the decomposition of m in the case of even n. 

We see, that in all cases m(x) is a quadratic form and 
in all cases its quadratic part is the same and equals to 
the symmetric form S2 (fTTj) . Below the notion for equiva- 
lence of Bell inequalities will be introduced. The maximal 
quantum violation vj of equivalent Bell inequalities is the 
same. It was shown in Q that all inequalities that max- 
imize Vf are equivalent. From the considerations below 
we can conclude: a boolean function f S F n maximize Vf 
if and only if it is of the form 

n 

/(x) =a(x)©s 2 (x) = Co ®0W®0^j. (H3) 

i— 1 i<j 

Affine functions correspond to the trivial Bell inequalities 
(|79p which cannot be violated in quantum case. Adding 
S2(x) substantially changes their properties: quantum 
violations become the largest. 

I end this subsection by showing that the boolean func- 
tions (|113p . which correspond to Mcrmin inequalities, 
minimize the uncertainty relation (j6~4")) . i.e. U(f) = 1 
for all functions of the form (|113[) . Since the uncertainty 
is invariant under equivalence of Bell inequalities, it is 
sufficient to check the equality U(f) = 1 only for one 
boolean function / of the form pi3p . Let us take s 2 . It 
is easy to see that 



NW S2 = 



if n is odd, 
if n is even. 



(114) 



Let us calculate NA S2 . We have 



(115) 



xev„ 



where £L- 



e 



We see that A S2 ^ only if u = 0, 
from which it follows that Ui = Ui = const. So 

there are at most two possibilities: u = (0, . . . , 0) and 
u = (1, . .., 1). For an odd n the number n — 1 is even 
and both possibilities give A S2 (u) = 0, from which we 
have NA S2 — 2, while for an even n only the former one, 
u = (0, . . . , 0), leads to A S2 (u) = 0, and in this case we 
have NA S2 = 1. In both cases we have U(s2) = L 



IV. CLASSIFICATION OF BELL 
INEQUALITIES 

The number of Bell inequalities 2 2 grows extremely 
fast with n. On the other hand, many of them are very 
similar, for example, differ in sign, and due to this such 
inequalities have similar properties (for example, they 
have the same maximal quantum violation). It makes 
sense to introduce the physically motivated notion of the 



equivalence of Bell inequalities and study only the classes 
of incquivalent inequalities. 

Two Bell inequalities (and the two corresponding 
boolean functions) are said to be equivalent if they can 
be obtained from each other by applying the following 
three kinds of transformations (any number of times, in 
any order): 

(i) permuting subsystems, 

(ii) swapping the outcomes of any observable, 

(iii) swapping the observables at any site. 

Let us study these three kinds of transformations and the 
equivalence of Bell inequalities in more details. 



A. Transformation of the first kind 

The permutation of subsystems, corresponding to a 
permutation ir £ S n (S n is the symmetric group), in 
terms of Bell inequalities reads as: with a Bell inequal- 
ity (|78[) corresponding to / € F n we associate another 
inequality with coefficients Wf(u) = Wf(iru), where 
7ru = 7r(ui, . . . ,u n ) = (u w -i(i),...,u w -i(„)). It is easy 
to see that the function W/(7ru) satisfies the condition 
(f24|) and due to this the function /' exits. Using the in- 
verse Walsch-Hadamard transform (f2"2"|) we can find the 
relation between / and /': 



(-1) 



/'(*) 



1 W)(7TU) 



uev„ 



(-l) { ^ u) W f (nn) 

= 4 E(- 1 ) <mc,u>M/ /( u ) = (- 1 ) / 



(116) 



(ttx) 



Here we used the fact that (7rx, ttu) = (x, u) for all x, y € 
V n , it € S n and that 



E = E w 

uev„ uev„ 



(117) 



We see that a permutation of subsystems is a permuta- 
tion of the arguments of the corresponding boolean func- 
tion: /'(x) = /(ttx). 

Let us define the map p n : T n — » T n via 



(P ff /)(X)=/(7TX). 



(118) 



Then any transformation of the type ((TJ) is p v for the 
appropriate tt. Since 7r(crx) = (<T7r)x, we have 



(P*P<rf)(*) = (Po-/)(7I" x ) = /(o-(7I" x )) 
= /((tTO-) x ) = (P7r<r/)(x), 

from which we get the equality 

PttPct — Pna ) 



(119) 



(120) 
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which is valid for all tt, a G S n . In other words, the map 
7r — > is a representation of S„ in i 7 ^. 

In terms of the states \ Wf) the transformation of this 
kind is given by the operator P^, defined via: 



PJx) = Itt^x) 



(121) 



It is easy to see that W/' (u) = W/ (7ru) if and only if 

\W r )=P v \W f ). (122) 

This operator is nonlocal — it permutes subsystems. 
Note that P v P a = P nrr . 

B. Transformation of the second kind 

Swapping the outcomes of the observable A(0) results 
in the following relations on the coefficients of Bell in- 
equalities: 



(123) 



U/'S. U/l. ...()....;. 

li/i.. U/;....i....;. 

Similarly, swapping the outcomes of -4,(1) we have 

w r (...,o,.-.) = w f (...,o,-.-), ( 

W fl (...,l,...) = -W f (.. .,i,...)- 1 ' 

These relations can be written as 

Wf.(u)=±(-iy*W f (u), (125) 

where + corresponds to Aj(l) and — to Ai(0). Swapping 
the outcomes of several observables Ai(u) with indices 
i € I C A/"„ = {1, . . . , n} has the following effect: 



W f ,(u) = ±(-l)^^W f (u) 



(126) 



where all components of e/ € y„ are zero except those 
with indices in /. Any boolean vector y <E V n can be rep- 
resented in the form e/: for any y € V n there is / C M n 
such that y = ej. Due to this we characterize transfor- 
mations of the second kind by a boolean vector yey„, 
and write the relation Q126P as 



W fl {u) = ±{-l) {y ' u) W f (u). 



(127) 



From this it is easy to find the relation between / and 
/' explicitly. For the sign + in (|127p we have 

(-!) /,(x) 4E (-l) (xey ' u) W?(u) = (-1)>C*»0, 

(128) 



so that 

/'(x)=/(x©y). (129) 
Analogously, for the sign — in Q127P we have 

/'(x) = /(x©y)©l. (130) 



Let us introduce the map 8 : F n — > F n and for any 
y G V n the map s y : F n —> F n via 

(i/)(x) = /(x)el, M(x) = /(xffly). (131) 

From the discussion above it follows that any transfor- 
mation ([u]) is s y or 6s y for an appropriate y £ V n . 

Consider the single qubit quantum gate with the ma- 
trix 



Z 



1 



(132) 



Due to the equality 

(u\zy\W f ) = (-1)<*»>W)(u) (133) 

the relations (|129[> and (|130|) are equivalent to the fol- 
lowing ones: 

\Wf,)=±Zy\W f ). (134) 
This transformation is local. 

C. Transformation of the third kind 

Swapping the observables of the z-th site, A;(0) *-> 
Aj(l), is expressed as 



(135) 



w>(... I o,...) = w>(...,i,...). 

W>/(... I 1,...) = W>(...,0,...), 

or in a more compact form as 

W7,(u) = W/(u©e 4 ). (136) 

Swapping the observables on several sites with indices 
J Q Nn reads as 



W f >(u) = Wf{u®ej). 



(137) 



\/(x)©(x,z} 



(139) 



As before, we characterize the transformations under 
study by a boolean vector z = e,/ and write the rela- 
tion (fT57| as 

W>(u) = W)(u®z). (138) 

It is easy to find the relation between / and /' explic- 
itly: 

(-!) /,(x) = 4 £(-l) (x ' u> W7(u©z) 

uGV„ 

4e (-i) (x ' uez ^/(u) = (-i) j 

u£y„ 

For any z € V„ let us introduce the map t z : F„ — > F n 
via 

(i z /)(x) = /(x)0(x,z). (140) 

We see that any transformation (jiu|l is t z for an appro- 
priate z £ V n . 

The relation (|138jl is equivalent to the following one: 

\W r )=X*\Wf). (141) 

This is also a local transform. 
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D. Relations between the transformations 

Let us analyze relations between the maps p^, 5, s y and 
t z . It is easy to see that S commutes with all the other 
maps and satisfies the relation 5 2 = 1. Composition of 
two s y or two t z is again a map of the corresponding 
kind: 



^y^y' ^y©y' 3 t z t z < t z 



(142) 



We see that the maps s y form a group, isomorphic to the 
additive group V n ; the same is true for the maps t z . 

Straightforward calculations show that the following 
relations are valid for all y,z e V„ and tt € S n : 

PirSy = S n -l y p„, Pirt z = t^-l z p n (143) 

Let us prove only the second relation. We have 

(p*i,/)(x) = (iz/)(7rx) = /(ttx) (ttx,z). (144) 
On the other hand, we have 

(^-^x/)(x) = /(ttx) © (x, ir- 1 *). (145) 
Taking into account that 

(7TX, z) = (n~ 1 TTX, 7T _1 Z) = (x,7T _ z), (146) 



we see that the last equality in (|143|) is valid. 

The maps s y and t z commute modulo the map 5: 



s y t z — 



t z s y if (y,z) = 0, 
St z s y if (y, z) = 1. 



In fact, we have 



(Vz/)(x) = (*z/)(x®y) 



/(x0y) © (x,z) © (y,z) 



(147) 



(148) 



Since 



/(xffly)ffl(x,z) = (t z s y /)(x), (149) 
the equality (|147|) is valid. 

E. The group Q n 

All the maps p n , 6, s y and t z are invertible and they 
generate a subgroup Q n of the symmetric group Sp n (the 
group of all one-to-one transformations of F n ). From the 
relations between p^, 5, s y and t z it follows that any 
element a G Q n can be represented as 



a = et z p^Sy 



(150) 



in an unique way, where e is either id : F n — * F n or 5. In 
particular, \Q n \ ~ 2 2n+1 nl. The element a (|150[) acts on 
a function / G F n as 

(a/)(x) = { /( " X0y)0(X ' Z> if £ = id ' (151) 

l/(7rx©y)8(x,z)©l if e = 5. 



This action can be written in the following form: 

(a/)(x)=/(7rxey)ffio(x), (152) 

where a £ A n is an affine function. Now we can for- 
mulate the equivalence of boolean functions (and the 
corresponding Bell inequalities) as follows: two boolean 
functions /, /' G F n are equivalent if there is a G Q n 
such that f = af. It is easy to see that it is really 
an equivalence on the set F n . Note that / ~ /' if 
and only if / © S2 ~ /' © S2 from which it follows that 
S : F n /~-^ i*n/~j [/] l— ► [/©S2] is a nontrivial involution 
on the set of the equivalence classes F n /~. The equiva- 
lence class [0] is the class of trivial inequalities and the 
class S([0]) = [52] is the class of the Mermin inequalities. 
The group Q n can be identified with the set 



Q n = {id, 8} X V n x S n x V n , 
equipped with the product: 

(e,z,7i",y)(e',z',7i"',y') 

= {e ee', z © 7r - V, tttt', Tr'y © y'), 



where £q is defined as 



id if (y,z')=0, 
8 if <y,z') = l. 



The unit element e of this group is 
e = (id, 0, id, 0) 



(153) 



(154) 



(155) 



(156) 



where the first id : F n — > is the identity map on i 7 ^ 
and the second id : Af n — > A/"„ is that on M n . The inverse 
(e, z, 7r, y) _1 reads as 

(EjZ.tt^)" 1 = (££,7rz,7r" 1 ,7r" 1 y), (157) 
where e is defined via 



id if (y, 7rz) = 0, 
6 if (y,7rz) = 1. 



(158) 



The elements (e, z,id, 0) form a normal subgroup 

U n < The elements (id, 0, 7r,y) form a (not nor- 
mal) subgroup J'n, usually referred to as Jevons group. 
The group Q n is the semidirect product ofU n 1%> and J n : 

Q n =U n V\J n . (159) 

The elements (id, 0, id, y) form a normal subgroup of J^ n , 
which is isomorphic to product of n copies of the cyclic 
group C2 of second order. The elements (id, 0, w, 0) form 
a (not normal) subgroup of J n which is isomorphic to 
the permutation group S n . The Jevons group J n is the 
semidirect product of C% and S n : 



Un 



C2 ^\ S n 



(160) 
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so we have the following decomposition of Q n : 

g n =u^\{c%\s n ). (i6i) 

In terms of vectors \Wf) the equivalence can be ex- 
pressed as: two boolean functions /, /' £ T n are equiv- 
alent if and only if there are 7r £ S n and y, z £ V n such 
that 



\W f .) = ±P ir Z*X'\W f ). 



(162) 



This gives another interpretation of the action of the 
group Q n on the set F n . 

The classes of equivalence 
of Bell inequalities for small 
number of qubits are shown 
in Fig. [5] For n — 2 there 
are N2 = 2 classes of equiva- 
lence, N3 = 5 and N4 = 39. 
The small squares, outlined 
in this figure, are shown in 
Fig. [6] They have the struc- 
ture similar with that of the 




Fig. 5(b) 



The blue subsquares (which form a circle) 
correspond to Mermin inequalities, and the orange ones 
(which form a cross) correspond to trivial inequalities. 
The same picture is repeated in these subsquares, show- 
ing a kind of fractal behavior. The Mathematica code 
with which these figures were obtained is presented in 
the Appendix A. 

The figures [5j as well as the figures [3] and |4] were ob- 
tained numerically. The analytical expressions for the 
number of equivalence classes and for the maximal quan- 
tum violation Vf of the equivalence class defined by a 
boolean function / £ F n are unknown. In the next two 
subsections a partial approach to this problem will be 
presented. 



F. Polya theory 

Consider two finite sets A and Y. The notation Y x 
is used for the set of all the functions / : A — > Y. The 
set Y x is finite and has |Y"|' X ' elements, which motivates 
the notation. Let G be a finite group and let a : G — ► Sx 
be a homomorphism of G to the symmetric group on X, 
i.e. an action of G on X. This means that for any g £ G 
the map 



a g = a(g) : X -> X 
is a permutation of X and the relation 



(163) 



l gg 



(164) 



is valid for all g, g' £ G. 

Two functions f,f f :X—*Y are said to be equivalent 
with respect to the action a, f ~ Q /', if there is g £ G 
such that /' = / o a g , i.e. if the relation 

f'{x) = f{a g x) 



is valid for all £ £ A. What is the number N of equiv- 
alence classes? For example, if the action a is trivial, 
i.e. if a g = idx for all g £ G, then / ^ a f if and only 
if / = /', which means that each equivalence class con- 
sists of a single element and there are IFI^ equivalence 
classes. 

To calculate the number of equivalence classes we need 
the notion of the cycle index of the group G. Each el- 
ement g £ G defines an equivalence ^ g on the set X 
via 

x ~ x' if x' = a g {x). (166) 
Each equivalence class has the form of a cycle 



z g (x) = {x,a g (x),...,a k(x) 1 (x)} 



(167) 



of some length k(x) (so that a g ^ (x) — x). Different 
cycles z g (x) do not intersect. For any fixed g £ G the set 
X can be represented as a disjoint union of cycles 

X = {Jz g (xi), z g (xi) H z g {xj) = if i^j. (168) 



Let Ck(g) be the number of the cycles of the length k in 
this decomposition, fc=l,...,n=|A|. It is clear that 
the numbers C\{g), . . . , c n (g) satisfy the following relation 



k=l 



kc k {g) 



(169) 



The cycle index of G with respect to the action a is the 
polynomial Zq(x\, . . . , x n ) of n variables defined as 



Z G (xi,...,x n ) = -n^X^i 1 
11 gee 



(9) ~.c»(g) 

■ ■ ■ Jj n 



(170) 



A very special case of the Polya theorem is the follow- 
ing statement: for the number of the equivalence classes 
( with respect to the equivalence (|165p ) we have 



N = Z G (\Y\,...,\Y\ 



(171) 



Now consider a more general situation, where there is 
an action not only on X but also an action f3 : H — > Sy 
of the group H on Y. The equivalence of functions given 
by (|165p can be extended as follows: f,f':X—*Y are 
said to be equivalent, / ~ /', if there are g £ G and 
h £ H such that the diagram 



X -^H-> A 

V \' 

Y Ph ) Y 
is commutative, i.e. if the relation 

(165) P h f'{x) = f(a g x) 



(172) 



(173) 



(a)n = 2 



(b)n = 3 




(c)n = 4 



FIG. 5: The classes of equivalence of Bell inequalities for small number of qubits. 
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FIG. 6: The subsquares. 



is valid for all x £ X. In this case the number N of 
equivalence classes can be expressed in terms of the cycle 
indices Zg(x\, . . . , x n ) and Zn{yi, • ■ ■ , y m ) as follows: 

N = Z G (^-,...,^-)z H (e s ^,...,e ms ^), (174) 

where Si(z) = Y^j=i z j ana - the derivatives are taken at 
Z\ = . . . = z n = 0. This expression can be transformed to 
another form, sometimes more suitable for calculations: 

N =W\J2 Z G(d(h), ci(ft) + 2c 2 (h), . . .), (175) 

' ' h£H 

where the i-th argument is equal to 

»-- ,:/ ' K ( 176 ) 

J I > 

It easy to see that the previous case (where there is a 
action only on X) is a special case of this more general 
situation when the action of H on Y is trivial. In fact, 
in such a case we have 

Cl {h) = \Y\, c k (h)=Q if k>l, (177) 



for all h e H, and from ([1751) we get the relation ([T7T]l : 



N = W\ ^ z{in ■ ■ ■ ' = z(in ■ ■ ■ ' |r|) - (178) 



If the equivalence on the set Y can be represented 
in the form (|173p then the number of equivalence classes 
can be calculated according to (|175p (or (|17ip ) . Unfortu- 
nately, not any equivalence on Y x can be represented in 
the form (|173(1 . In such cases one must find other ways 
to solve the problem. 

G. Classification with respect Jevon's group J n 



In our case X = V n and Y = Z2. The equivalence 
(|152p cannot be represented directly in the form (|173p . 
so let us start with a simpler case. Consider the action 
of J n on V n , which is the reduction of the action of Q n 
to its subgroup Jn ■ The cycle index of the group J n 
was calculated in [HI and it is given by the following 
complicated expression: 



Z Jn (x 1 ,...,x 2 n) = Y^ 



J! c,!(2i) c * ^=l 



■® n 



n 

d\1i,d\i 



) 



(179) 
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n 


Zj n (xi, . . . ,xa») 


N n 


N n 


I 


1 / 2 . \ 

2 W + X2 i 


2 


I 


2 


+ 3^2 + 2x^X2 + 2x 4 ) 

8 


4 


2 


3 


i(a;f + 13^2 + 8x1X3 + 8x 2 x 6 
48 

+ 6x1X2 + 12x4) 


14 


5 


4 


1 / IS , tl 8 , JO 2 3 

— ■ (xi + 51x 2 + 48xix 2 x 4 
384 

+ 48xg + 12xf A + 84x1 

+ \2x\x\ + 32x*X3 + 96x2X5) 


222 


39 


5 


-^_(xf + 231x 2 6 + 20x; 6 x« 
+ 520x| + 80x?X3 + 720x2Xg 

+ I6OX1XIX3X6 + 320x 4 ! x?2 

+ 2<XSx\x\x\ + 480x| + 24OX2X4 

+ 60xfxa 2 + 384xiXg + 384x 2 x? ) 


616126 


22442 



TABLE III: Cycle index Zg n for small n. 



where gcd(p, q) and lcm(p, q) are the greatest common 
divisor and the least common multiple of integers p and 
q respectively. For two monomials the cross-product is 
defined via 

k I 

( X ;i ...x^)® {xz 1 ...o = n n Kt ® x ?/)> ( i85 ) 

»=i i=i 

and then extended for arbitrary polynomials by bilin- 
carity. For example, let us calculate the cross-product 

{xlxD ® {xix\): 

{x\x%) <8> (xjxl) = (x\ <E> xf)(x\ <g> x\){x\ <g> x\)(x\ <g> x\) 



„2-4-l„2-5-l„3-4-l„3-5-2 



„8 40 12 
X3X4 Xg . 



(186) 



The cycle index Zj n for small n is shown in the second 
column of the tabic IIII1 

Now let us add an action on the set Y = Z2. Let 
H = C2 = {1,t} be the cyclic group of the second order. 
The action (3 we define as: /3\ = idz, and (3 T being the 
logical NOT, (3 T (0) = 1 and /3 T (1) = 0. For cycle lengthes 
we have 



ci(l) 



c 2 (l)=0, Cl (r) = 0, c 2 (r) = l. (187) 



According to (|175[) for the number of equivalence classes 
we have the following expression (see also fill]): 



where the sum is over all vectors c = (ci, . . . , c„) with 
nonnegative integer components such that 



(180) 



k=l 



The function a(/c) is defined for any integer k ^ 1 via 



d\k 



dj' 



(181) 



and b(2k) is defined for positive even integer argument 
2k, k ^ 1, via 



x»)-5 s 2i,2 "(¥)- 

with fJ,(m) being the Mobius function 

1 if m = 1, 

MM = <( (-l) fc if m=p 1 ...p k , 
in other cases, 



(182) 



(183) 



where pi, . . . ,pk are different prime numbers. 

The cross-product is defined as follows. For powers of 
variables we have 



^ m _ nmgcd(p,g) 
x p w x q — x lcm ( p ^ , 



(184) 



Z G (2,...,2) + Z G (0,2,0,2,...) 



(188) 



These numbers are shown in the third column of the ta- 
ble Mil The last column of the table shows the num- 
ber N n of equivalence classes with respect to the equiv- 
alence under study. The number N§ was taken from 
www.ii.uib.no/ larscd/boolcan/. The numbers N n for 
n > 5 arc unknown. 



V. CONCLUSION 

In conclusion, the relation between the boolean func- 
tions theory and the general Bell inequalities for n-qubits 
is established. The classification of Bell inequalities with 
respect to the Jevons group is obtained, which is a weaker 
result then the problem posed in Q. Nevertheless, to 
my knowledge it is the only approach to the more gen- 
eral classification. This approach is based on the works 
[13, EH done for computer logic circuits theory, which 
shows the connection between quite different problems 
- qubit system description and computer logic circuit 
design. 

There are still many unsolved problems. Two the most 
important ones are: 

(i) classification of Bell inequalities (or boolean func- 
tions) with respect to the group Q n , 
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(ii) characterization of the maximal quantum violation 
v f of a given boolean function / in terms of proper- 
ties of /, in particular, finding the relation between 
vj and the nonlinearity Nf. 

And, of course, a very interesting question — how the 
ideas from different applications of the boolean functions 
theory (not only to cryptography or computer logic cir- 
cuits design) can be used in quantum information theory. 

APPENDIX A: THE CODE 

In this Appendix I present the Mathcmatica code Figs. 
[2 H H and Table [En] were obtained with. First of all, 
the set V n can be coded as 

in[i]:= V[n_] := Tuples [{0, l},n] (Al) 



where J is the list of all the multi-indices (ii, . . . , ik) for 
which the coefficients c^,...^,. 7^ and c the list of these 
coefficients. Note that the Length-es of c and J must be 
the same. For example, the expression 

f [c, Subset s [Range [n], {m}]] (A3) 



Similarly, the expression 

f [c, Subsets[Range[n],m]] (A5) 



This method of calculation the Walsch-Hadamard trans- 
form is very simple and quite ineffective. In Appendix B a 
much better method is presented. The only disadvantage 
of that method is the fact that it works much faster only 
when it is necessary to calculate all the numbers Wf(u), 



The Mathcmatica function IntegcrDigits[n, b, 1] gives the 
list of the base-6 digits of n, padding it on the left if 
necessary to give a list of length I. Using this function, 
one can code V„ in another way as 



V[n_] := IntegerDigits[#,2,n]&/@Range[0,2 n - 1] 



This code works a few times more slowly, but this ap- 
proach can be useful if it necessary to construct only 
some part of V n , not the whole set V n . We know that 
F n ~ V2™, but since 2 2 > 10 19 , the simple definition 
F[n_] := V[2 n ] will not work for n ^ 6. 

A general boolean function ([8]) can be coded as follows: 



(A2) 



gives a homogeneous polynomial of the degree to (except 
the case when all elements of c are zero; when all ele- 
ments of the list c are 1 then (|A3|) is the m-th symmetric 
polynomial (flT|) ; 



(A4) 



is a polynomial of the degree not greater then to. 
The Walsch-Hadamard transform can be coded as 



(A6) 



u G V n simultaneously, and it cannot be applied to calcu- 
late only one number Wf (u) for a given u € V n . WHT is a 
functional: WHT[f ] is a function and WHT[f ] [{ui, . . . , u n }] is 
its value W/(u) at u = (ui, . . . , u n ). The inverse Walsch- 
Hadamard transform can be coded as 



I 

in[2].-= f[c_,S_] := Function[x,Mod[c.(Apply[Times, #]&/@ (Part[x, #]&/@ S)),2]] 

I 



I 

f [Table [1, {Binomial[n, m]}], Subsets [Range [n], {m}]][{xi, . . . , x n }] = s m (xi, . . . , x n ) 

I 



I 

in[3]:= WHT := Functional, Plus @@ (Function[x, (-l) #[x]+(x - u) ]/@ V[Length[u]])]& 

I 



/n[4] : = WHIT := Function[x, (1 - 2" Length[x] Plus @@ (Function[u, (-l) x ' u #[u]]/@ V[Length[x]]))/2]& (A7) 
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Like WHT, it is also a functional. be coded as 

The autocorrelation (|46|) of two boolean functions can 

I 

/„[5]:= A := Function[u,Plus @@ (Function[x, (_ 1 )#iM+#2[Mod[#,2]&/<a («+u)]]/@ V [Length [u]])]& (A8) 

I 

It is a functional of two arguments: the expres- (ui,...,u n ). 

sion A[f , g][{ui, . . . ,u n }] is the value A/ jS (u) at u = The uncertainty (|65[) can be coded as 



in[6]:= U[n_] := Function[f,2~ n Count[WHT[f]/@ V[n], _?(#! = 0&)]Count[A[f, f]/@ V[n], _?(#! = 0&)]] (A9) 

I 

For a function / £ F n the expression U[n][f] gives the we need a way to define / given its number B n (f), ^ 
uncertainty U(f) of /. To calculate U(f) for all / £ F n B n (f) < 2 2 . The following code solves this problem: 



ln[7]~ itof[n_, B_] := Part[lntegerDigits[B, 2, 2 n ], 2 n - FromDigits[#, 2]]& (A10) 

I 

Given B < 2 2 , the expression itof[n, B] is the ize boolean functions with respect to their uncertainty 
corresponding boolean function, to which one can ap- using the code 
ply the syntaxis itof [n, B][{x 1 , . . . ,x n }]. We can visual- 



in[8]~ ut[n_] := Partition[Table[U[n][itof [n,B]], {B, 0, 2 2 - 1}], 2 2 ] (All) 

I 

In Mathematics this table can be immediately plotted For n = 2 it is 1 (all 16 boolean functions have the same 

with the ArrayPlot function, what was doe for the case uncertainty 1), for n — 3 it is 2 (the values are 1 and 

of n = 4, for the other two cases I used a simple script to 8) and for n = 4 it is 4 (the values are 1, 35/8, 8 and 

generate pstricks code from this table and then compiled 16). Unfortunately, these are all the values of n for which 

it with WF$i.(pstricks code produced huge pictures in this simple code works, for larger n another technique is 

the case of n = 4). In this way Fig. [3] was obtained. The needed, 
code 

ut [n] / /Flatten/ /Sort / /Split / /Length (A12) 

Now I will show how the Fig. [2] was obtained. The key 

gives the number of different values of the uncertainty. point is the function 



ln[9]~ d[n_, m_] :— Function[c, FromDigits[f [c, Sub sets [Range [n], m]]/@ V[n], 2]] 
/@ V[Plus @@ (Binomial [n,#]&/@ Range [0, m])] 



I 

which returns the numbers B n (f) of the boolean func- tions / £ F n of degree ^ m. To visualize boolean func- 
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tions with respect to their degree let us create a list The table Partition[dt, 2 2 " ] gives the desired visual- 

dt = Table[n, {2 2 }];. Then we can fill it in as ization. 

For[i = 1, i ^ n, i + +, dt [[d[n, n — i] + 1]] = n — i] Now let us discuss the equivalence of Bell inequalities. 

(A14) The map p v (|118p can be coded as 



in[w]-.= p[pi_] := Function[x, ^[Permute[x, InversePermutation[pi]]]]& 

I 



(A15) 



The expression p[pi] is a functional: p[pi][f] gives the (p ff /)(x) at x = (x\, . . . ,x n ). The maps 5, s y (|131[) and 
function p-nf, and p[pi][f][{xi, . . . ,x n }] gives its value t z (|140[) can be coded as 



In[U] 
In[12] 
/n[13] 



= 8 := Function[x, 1 - #[x]]& 

= s[y_] := Function[x, #[Function[z,Mod[z, 2]]/@ (x + y)]]& 
= t[z_] := Function[x,Mod[#[x] + x.z, 2]]& 



(A16) 



r 



Let us illustrate the relations between the maps under 
study. As an example, consider the first relation (|143| . I 
show that this relation is valid by applying the maps from 
both sides to all boolean functions / G F„ and comparing 



J 



the results. To do it, we need operations which produce 
the lists of values {(p w s y /)(x)} and {(s y p,r/)(x)}, x G 
V n given the list {/(x)}, x G V n of values of / G F n . The 
code 



/n[i4] : = ps[n_,pi_,y_] 
Jn[i5]:= sp[n_,pi_,y_] 



Composition[p[pi], s[y]][Function[x, Part[#, FromDigits[x, 2] + l]]]/@ V[n]& 
Composition[s[y],p[pi]][Function[x,Part[#,FromDigits[x, 2] + l]]]/@ V[n]& 



(A17) 



is a solution to this problem. The expression 



J„[16]:= ps[3, {2, 3, 1}, {1, 1, 0}]/@ V[2 3 ] - sp[3, {2, 3, 1}, {1, 0, l}]/@ V[2 3 ]//Short 

I 

produces a list of zero-lists, as expected. This example Given a number ^ B < 2 2 the function 

clearly demonstrates that the first relation (|143| is valid 
(in the case of n = 3). 



(A18) 



{{0,0,0,0,0,0,0,0},<<254 >>, 
{0,0,0,0,0,0,0,0}} 



(A19) 
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/n[i7]~ e[n_, B_] := Replace[FromDigits[Map[#, V[n]] , 2]& 

/@ Function[f , Apply [Function[{e, y,pi, z}, Composition^, t[z],p[pi], s[y]] 
[Function[x, Part [lntegerDigits[f , 2, 2 n ], FromDigits[x, 2] + 1]]]], #]& (A20) 
/@ Tuples[{{Identity, <5}, V[n], Permutations [Range [n]], V[n]}]][B] 
//Sort//Split,{(i.)..} i,l] 



produces the list of numbers which correspond to the 
functions, equivalent to the one corresponding to B. 
With this functions it is easy to get the equivalence 
classes. Let us illustrate the general idea by the case 
of n = 3. We start with 5 = 0: e[3, 0]//Short produces 



{0, 15, 51, << 10 >>, 204, 240, 255} 



(A21) 



It is the first class of equivalence (which contains 16 el- 
ements). The smallest number which is not in this class 
is 1; e[3, l]//Short produces 



{1,2,4, << 122 >>, 251, 253, 254} 



(A22) 



It is the second class of equivalence (which contains 128 
elements). The smallest number which is not in cither 
class found so far is 3; e[3, 3] //Short produces 



{3, 5, 10, << 42 >>, 245, 250, 252} 



(A23) 



It is the third class of equivalence (which contains 48 
elements). The ext number to try is 6: e[3, 6] //Short 
produces 



{6, 9, 18, << 42 >>, 237, 246, 249} 



(A24) 



It is the fourth class of equivalence (which also contains 
48 elements). The next number is 23; e [3, 23] //Short 
produces 



{23, 24, 36, << 10 >>, 219, 231, 232} 



(A25) 



it is the fifth class of equivalence (which contains 16 ele- 
ments). Since 16+128 + 48 + 48 + 16 = 256, we exhausted 
all numbers (boolean functions) which means that in the 
case of n = 3 there are 5 classes of equivalence, they are 
sho wn in Fig. |5(b)| The same approach allows one to get 



Fig. 5(c 



The maximal quantum violations can be found as fol- 
lows. Let us introduce the functions 

In[i8]~ S[ip.] := {Cos[(p],I S±n[tp]} 
/n[i9]:= S^.List] := Apply[Times, #]& (A26) 
/@ Tuples [S/@ <p] 

The maximal violation (1861) can now be coded as 



(A27) 



/n[20]~ v[f _,¥?_] := Abs[((-l) f[#1 & 

/@ V[Lengt%]]).S[y>]] 
in[2i]-.= v2[B_, ip_] := v[itof [Length[^], B], 



The quantity v2[B, </?] give the maximal quantum vio- 
lation Vf of the function / corresponding to the num- 
ber B. Let us illustrate the calculation of maximal 
quantum violation for Mermin inequalities. The Mer- 
min inequalities (whose coefficients are given by (|9ip and 
(|52l ) can be coded as (using the standard package Alge- 
bra 'SymmetricPolynomials ) 

m[x_] := Mod[SymmetricPolynomial[x, 2], 2] (A28) 

The maximal quantum violation of Mermin inequalities 
can be found as 



(A29) 



/n[22]:= FindMaximum[v [m, {(p±, . . . , fn}], 
{<p u </??},... ,{<Pn,<pl}] 

For example, the code 

FindMaximum[v[m, {<p u <p 2 }],{tpi, 1}, {f2, 1}] (A30) 
produces 

{1.41421, {ifi -> 0.785398,v9 2 ^ 0.785398}} (A31) 
the code 

FindMaximum[v[m, {tp±, tp 2 , ^3}] , 

{tp±,l},{<p 2 , 1},{^3, 1}] 

produces 

{2.,{ip t —> 0.785398,(^ 2 ^ 0.785398, tp 3 - 



(A32) 



■•}} 

(A33) 

in full agreement with the relation (|87|) . Since Vf is the 
same for equivalent functions, using the classes of equiv- 
alence calculated before, it is easy to get Fig. [4] 

Now I will show how the table Mil was obtained. To 
calculate the cycle index Zj n (x±, . . . , x 2 rC) we need to cal- 
culate the sum ()179p . The vectors c for a given n can 
be obtained using the standard package Combinatorica, 
containing the function Part it ions [n], which returns the 
list of partitions of n. Any partition p of this list has the 
following form 



p = {m, . . . ,ni,n 2 , . . . ,n 2 ,n 3 , . . .}, 



with ni > "2 > ^3 > 
relation 



and kxTii + ^2^2 



,h,.. .} 



(A34) 



n. The 



(A35) 
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gives a one-to-one correspondence between the partitions other components are zero. This correspondence can be 
of n and the vectors c. Here ki is on n^-th place and the realized in Mathcmatica via 



/n[23]:= ct[p_] := Module[{c = Table[0, {Plus @@ p}]}, 

Set[Part[c,#[[l]]],Length[#]]&/@ Split |p]; 

r i (A36) 
Return[c] 



I 

Then we need to define the cross-product. We introduce the p-th independent variable. The code 
the object var[p,n] which represents the n-th power of 



/n[24]~ cp[var [p_, n_] , var [q_, m_]] := var[LCM[p, q], n m GCD[p, q]] (A37) 

I 

reproduces the definition (|184[) . To extend the cross- definitions (the order in which they are given is impor- 
product for arbitrary polynomials we need the following tant): 



in[2b]-.= cp[v_, c_] := c v/; NumericQ[c] 

7n[26]:= cp[c_, v_] := c v/; NumericQ[c] 

/n[27]:= cp[vl_ + v2_, v3_] := cp[vl, v3] + cp[v2, v3] 

/n[28]:= cp[v3_, vl_ + v2_] := cp[v3, vl] + cp[v3, v2] 

Jn[29]:= cp[v3_, vl_v2_] := cp[v3,vl] cp[v3,v2] 

/n[30]:= cp[vl_v2_, v3_] := cp[vl,v3] cp[v2,v3] 

Jn[3i]:= cp[c_vl_, v2_] := c cp[vl,v2]/; NumericQ[c] 

Jn[32]:= cp[vl_, c_v2_] := c cp[vl,v2]/; NumericQ[c] 

I 



(A38) 



The following code reproduces the identities (x") 



J 



and 



n n+rn 



respectively: 



Jn[33]:= var/: Power [var [p_, n_], m_] 
in[34]:= var/ : var[p_, n_] var[p_, m_] 

I 



= var[p,nm] 
:= var[p,n + m] 



(A39) 



The code for the functions a (| 18 1(1 and b (|182p is obvious: 



/n[35]:= a[k_] := - plus @@ [2* MoebiusMu 
Zn[36].-= b[k_] := ^ plus ®@ ( 2 #/2 MoebiusMu 



k 

# 
k 

# 



&/@ Divisors[k] 
&/@ Complement [Divisors [k], Divisors[k/2] 



(A40) 
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A factor of of the sum (|179|) can be coded as 



/n[37]~ f actor[i_] := Times @@ (var[#, a[#]]&/@ Divisors[i]) + Times 
/@ Complement [Divisors [2i], Divisors[i]]) 



(var[#,b[#]]& 



r 



(A41) 



A term of the sum is coded as 



in[38]:= term[c_] := Module[{t}, 

t = Maplndexed[lf [#1 == 0, 1, If [#1 == 1, f actor[#2[[l]]], 

Fold[cp, f actor[#2[[l]]], Table[f actor[#2[[l]]], {#1 - 1 }]]]]&, c]; 
Foldfcp, First [t] , Rest [t]] 
]/(Times @@ (Factorial/® c)Times @@ Maplndexed[(2#2[[l]]) #1 &, c]) 



The whole sum (|179p is given by 



(A42) 



/n[39]:= cyclelndex[n_] := Plus @@ (term/® (ct/@ Part it ions [n])) 

I 



(A43) 



The number N n can be obtained as follows: 



in[40]-.= ub[n_] := -((cyclelndex[n] /. {var[p_, k_] — * 2 k }) 



(cyclelndex[n] /. {var[p_, k_] 
I 



If [EvenQ[p],2 k , 0]})) 



(A44) 



Then one can get the numbers form the ta- 
ble IIIII the code Table[ub[n], {n, 1, 5}] produces 
{2,4, 14,222,616126}. 



APPENDIX B: FAST WALSCH-HADAMARD 
TRANSFORM 



A much faster way of calculating the Walsch-Hadamar 
transform of a boolean function / G F n (i.e. calcula- 
tion of all 2™ numbers Wf(u), u £ V n ) is based on the 



following decomposition of the Hadamard matrix: 

n n 

H n = n ^-fc+i = n ® h ® ( bi ) 



k=l 



k=l 



which can be proved by induction. The relation y = M^x 
can be written as 



Uj2 k +i — x j2 k +i + x j2 k +2 k - 1 +i: 
Vj2 k +2 k - 1 +i — x j2 k +i — x j2 k +2 k - 1 +i> 



(B2) 



for i = 1, . . . , 2 fe " 1 and j = 0, . . . , 2 n ~ k - 1. The Walsh- 
Hadamard transform can be calculated according to (fTg|) 
as 



w f = M n . . .MxZf. 



(B3) 
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Starting with the vector x = Zf we calculate vectors 
Yk = Mk ■ ■ ■ M\Zf for k = 1, . . . ,n. Then y„ is w/. The 
algorithm in C is presented in the function whtl in the file 
whtl.c below. To turn whtl.c into a program usable from 
inside Mathematica another file, whtl.tm is necessary. 

whtl.tm 



: Begin: 

: Function: 

: Pattern: 

: Arguments : 

: Argument Types : 

: ReturnType : 

:End: 



whtl 
WHTL [x_List] 
{x} 

{IntegerList} 
Manual 



whtl. 1 



#include "mathlink.h" 
#include <stdlib.h> 
#include <string.h> 

// xlen is of the form 2™, returns n 
long logdong xlen) { 

long i = 1, m = xlen; 

while((m >>= 1) != 1) i++; 

return i; 

} 

void whtl (int *x, long xlen) { 

int *y = (int*)calloc(xlen, sizeof (int) ) ; 
long i, j, k, t, powl, pow2, indl, ind2; 
long n = log(xlen) ; 

for(k = 1; k <= n; k++) { 
powl = l«(k-l) ; 
pow2 = l«(n-k) ; 
for(j = 0; j < pow2; j++) { 



t = j«k; 

// in C arrays are numbered from 
for(i = 0; i < powl; i++) { 

indl = t + i; 

ind2 = indl + powl ; 

y[indl] = x[indl] + x[ind2]; 

y[ind2] = x[indl] - x[ind2]; 

} 

} 

if (k < n) 

memcpy(x, y, sizeof (int) *xlen) ; 



} 



MLPutlntegerList (stdlink, y, xlen); 
free(y) ; 



} 



int main(int argc, char* argv [] ) { 
return MLMain(argc, argv); 

} 

On a UNfX system the files are compiled with the fol- 
lowing command 



mcc — o WHTL whtl.tm whtl.c 



(B4) 



This command produces an executable file WHTL. To 
use it in Mathematica it must be installed as 

Install["/path/to/program/WHTL"]; (B5) 

Then it can be used as 

WHTL[{z 1; ...,z 2 n}] (B6) 



where all Zj, i = 1, . . . , 2™ are ±1. 
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